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ONSET  CRITERIA  FOR  STRUCTURE  IN 
MAGNETICALLY  CONFINED  PLASMA  EXPANSIONS 


I.  INTRODUCTION 

Recent  experiments  have  shown  the  rapid  development  of  structure 
during  the  early  phase  of  magnetically  confined  plasma  expansions 
[Bernhardt  et  al.f  1987;  Ripin,  1987].  Such  plasma  expansions  can  be 

t 

characterized  as  sub-Alfvenic,  i.e.,  <  V^a  where  is  the  initial 
expansion  velocity  and  V^a  is  the  Alfven  velocity  in  the  ambient  background 
plasma.  Under  this  condition  the  expanding  plasma  is  'stopped'  by  the 

r 

magnetic  field.  This  is  in  contrast  to  super-Alfvenic  expansions  where 
momentum  coupling  to  the  ambient  plasma  is  primarily  responsible  for 
'stopping'  the  expansion.  Another  feature  of  these  expansions  is  that  the 
directed  ion  Larmor  radius  (p.  «  ^q/2^  where  2.  is  the  ion  cyclotron 
frequency)  can  be  comparable  to  or  much  larger  than  the  scale  sizes  of 
interest.  For  instance,  in  the  cases  of  AMPTE  and  the  NRL  laser 


experiment,  p.  ~  Rg  »  Lr  where  Rg  is  the  maximum  radius  of  expansion  and 
Ln  is  the  scale  size  of  the  density  shell.  In  a  similar  vein,  one  should 
note  that  the  time  scale  («  *)  for  structure  development  in  these 
experiments  is  faster  than  an  ion  cyclotron  period  (2^-^),  i.e.,  »  >>  2^. 
What  makes  this  regime  interesting  is  that  it  cannot  be  described  by  ideal 
or  finite  Larmor  radius  (FLR)  MHD  theory  because  conventional  MHD  theory 
requires  «  <<  2^  and  p^  «  Ln< 

In  order  to  understand  the  gross  evolution  and  stability  properties  of 
magnetically  confined  plasma  expansions,  Hassam  and  Huba  (1987a, b)  recently 
developed  a  modified  set  of  one-fluid  MHD  equations.  These  equations  are 
valid  for  arbitrary  values  of  p^/Ln  (or  w/2j)  but  require  pg/Ln  <<  1  (or 
w/2e  «  1)  (i.e.,  the  electrons  are  strongly  magnetized).  Linear  stability 
analysis  indicates  that  a  fast  'unmagnetized  ion'  Rayleigh-Taylor 
instability  can  develop  during  the  early  phase  of  a  magnetically  confined 
plasma  expansion  when  p./Ln  »  1  (or  u  >>  2^).  The  instability  is  driven 
by  the  deceleration  of  the  plasma  shell  which  produces  an  effective 
gravitational  acceleration,  i.e.,  dV^/dt  ■  -  8eff?  acting  on  the  density 
gradient  associated  with  the  plasma  shell.  Moreover,  the  eigenmode 
structure  suggests  that  the  density  shell  should  exhibit  plasma  'clumping' 
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instead  of  plasma  'rippling'  (as  is  the  case  with  the  conventional 
magnetized  ion  Rayleigh-Taylor  instability).  Recent  2D  MHD  simulations 
also  support  this  type  of  behavior  (Huba  et  al.,  1987).  Both  of  these 
features,  rapid  development  and  density  clumping,  are  consistent  with 
observational  evidence  thereby  lending  support  to  the  Hassam  and  Huba 
model. 

The  purpose  of  this  paper  is  to  develop  a  rather  simple  set  of 
criteria  for  the  onset  of  structure  in  magnetically  confined  plasma 
expansions.  Ve  are  specifically  concerned  with  structure  generated  by  the 
unmagnetized  ion  Rayleigh-Taylor  instability  discussed  in  Hassam  and  Huba 
(1987a, b).  We  show  that  these  criteria  are  consistent  with  available 
experimental  data,  and  can  be  readily  tested  in  future  NRL  laser 
experiments. 

II.  THEORY 

A.  Magnetically  Confined  Expansions 

We  first  show  that  magnetically  confined  plasma  expansions  are  sub- 

/ 

Alfvenic.  We  estimate  the  'stopping'  radius  of  a  magnetically  confined 

expansion  by  equating  the  initial  kinetic  energy  of  the  plasma  with  the 

3 

magnetic  energy  in  a  volume  (4/3)nRB  where  Rg  is  defined  as  the  magnetic 
confinement  radius, 


1MV  2 

2  VdO 


t  &  V) 


where  Mn  is  the  mass  of  the  plasma, 


V^q  is  the  initial  debris  expansion 


velocity,  and  Bq  is  the  ambient  magnetic  field.  From  (1)  we  obtain 

/  2 

rb  *  3M0Vd0/B0 


We  compare  this  distance  with  the  equal  mass  radius 


which  is  defined  as  the  radius  of  the  sphere  which  contains  a  mass  of 
ambient  plasma  equal  to  the  mass  of  the  expanding  plasma.  For  the 
expanding  plasma  to  be  confined  magnetically  we  require  Rfi  <  R...  This 
leads  to  V^q  <  VAa  from  (2)  and  (3)  where  VAa  -  BQ/(4nnama)1/  is  the 


# 

Alfven  velocity  in  the  ambient  plasma.  Basically,  when  VJrt  «  VA  the 

dU  Aa 

expanding  plasma  is  stopped  before  it  sweeps  up  very  much  background 
plasma. 


B.  Magnetic  Deceleration 

The  driving  mechanism  of  the  unmagnetized  ion  Rayleigh-Taylor 
instability  is  the  deceleration  of  the  expanding  plasma  shell.  This 
deceleration  can  be  interpreted  as  an  effective  gravitational 
acceleration,  i.e.,  dV^/dt  -  -  ggff.  We  estimate  ggff  based  upon 
conservation  of  energy.  Ve  write 


?Movd2(t>  *  sr  rR3«>  -  iMovao2 


where  the  LHS  of  (4)  is  the  sum  of  the  kinetic  energy  of  the  expanding 
debris  and  the  swept-up  magnetic  energy  at  time  t  and  position  R,  and  the 
RHS  is  the  energy  at  t  *  0.  We  solve  (4)  for  V^t)  and  obtain 
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We  take  the  time  derivative  of  (5)  to  obtain 


dV,  Bn  7 

geff(t)  ’  '  dt-  =  2M^  R  (t) 


where  we  have  made  use  of  the  fact  that  «  dR/dt. 


We  also  note  that  the  position  of  the  debris  shell  can  be  determined 
as  a  function  of  time  from  (5).  Namely,  (5)  can  be  rewritten  as 
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Using  (7)  ve  can  determine  the  time  t  vhen  the  debris  shell  is  at  a  radius 
R. 

C.  Onset  Criteria  for  Structure 

The  unmagnetized  ion  Rayleigh-Taylor  instability  discussed  by  Hassam 

1/2 

and  Huba  (1987a, b)  can  only  be  excited  if  Yq  >  2^  where  Yq  *  (g  ^/L  ) 

That  is,  the  maximum  growth  rate  of  the  magnetized  ion  Rayleigh-Taylor 
instability  has  to  be  greater  than  the  ion  cyclotron  frequency  in  order 
for  the  ions  to  behave  as  unmagnetized  particles.  When  this  occurs  then 
the  unmagnetized  ion  interchange  mode  can  be  excited.  Making  use  of  this 
fact,  and  the  expression  for  derived  above,  we  then  arrive  at  the 

following  onset  criterion  for  structure  development  in  magnetically 
confined  expansions: 


or 


B  2  2 

0  R  t 1 ) 

2M0  Ln 


R(t)  >  Rc 


(8) 

(9) 


where  Z  is  the  charge  state  of  the  ions.  Note  that  R£,  the  critical 

radius  for  structure  development,  is  independent  of  the  magnetic  field. 

This  is  because  both  and  Yq  are  proportional  to  Bq.  We  do  not  expect 

that  structure  will  be  observed  at  R  =  R  but  that  this  is  the  radius  at 

c 

which  structure  can  begin  to  grow.  Since  it  will  take  a  number  of  e- 

folding  times  for  the  structure  to  be  observable,  we  anticipate  that 

structure  should  be  observed  at  R  ,  =  cxR  where  a  =  1.5  -  2.0. 

obs  c 

Another  constraint  on  structure  development  is  that  R0^s  =  «RC  <  RB* 
That  is,  the  magnetic  confinement  radius  must  be  large  enough  to 
accommodate  structure  development.  Alternatively,  if  RoBs  =  <*Rc  >  R0  then 
the  maximum  deceleration  of  the  plasma  shell  is  not  fast  enough  for  the 
ions  to  behave  as  unmagnetized.  Thus,  we  also  require 


(10) 


where  we  have  used  (2)  and  (9). 
experiment  we  can  rewrite  (10)  as 


With  regard  to  the  DNA/NRL  laser 


2 

vhere  V  is  the  laser  energy,  i.e.,  we  have  assumed  MqVjq  /2  =  W.  Thus, 
(11)  suggests  that  the  laser  energy  must  exceed  a  critical  energy  (Wc)  in 
order  to  observe  structuring  of  the  expanding  plasma. 

III.  APPLICATION 

We  now  apply  the  criterion  developed  in  Sec.  II  to  the  structure 
observed  in  AMPTE  and  the  DNA/NRL  laser  experiment.  We  rewrite  (9)  as 

4  Z  f  \^2 

R  >  Rc  =  1.3  x  104  [LnM0J  cm  (12) 

where  u  *  m^/nip,  m^  is  the  proton  mass,  and  we  have  evaluated  the  physical 

constants.  For  AMPTE  we  note  that  Z  -  1,  y  »  137,  and  Ln  e  4.0  x  10^  cm. 

We  take  Mq  *  N^ym^  where  Nq  is  the  total  number  of  ions  produced.  From 

Bernhardt  et  al.  (1987)  we  take  Nq  =  7.5  x  10^  ions  so  that  Mq  =  1.7  x 

10^  gm.  Thus,  for  AMPTE  we  find  that  Rc  =  80  km  and  anticipate  structure 

to  be  observed  at  Rol)S  -  «RC  =  120  -  160  km.  This  is  consistent  with 

observations  which  indicate  structure  is  evident  at  R  a  150  km.  For  the 

laser  experiment  we  take  Lr  =  0.5  cm,  Mq  =  2.0  x  10-7  gm,  y  =  28,  and  Z  = 

10.  We  choose  the  value  of  Z  based  on  the  NRL  HANEX  simulations 

(Mulbrandon  et  al.,  1987).  We  find  that  R  =  1.5  cm  so  that  structure 

c 

should  be  observed  at  R0jjS  =  2.2  -  3.0  cm.  This  is  consistent  with 

experimental  observations  for  both  Bq  =  1  kG  and  10  kG. 

We  can  also  rewrite  (11)  as 

11  3  2  Z^  (  \3/2 

W  >  Wc  =  3.3  x  10iA  «\l  ^  [LnM0J  erg  .  (13) 

For  the  laser  experiment  we  note  that  Wc  =  16  -  38  J  for  Bq  =  10  kG  and  W£ 
=  0.16  -  0.38  J  for  Bq  =  1  kG.  Since  W  =  30  J,  (13)  is  marginally 

satisfied  for  Bq  =  10  kG  and  easily  satisfied  for  BQ  =  1  kG.  We  would 
expect  that  the  unmagnetized  ion  Rayleigh-Taylor  instability  would  not 
develop  if  W  < 


We  now  present  a  series  of  plots  which  graphically  illustrate  the 

above  criteria.  We  consider  the  following  parameters  which  are  relevant 

to  the  laser  experiment:  Bq  =  10  kG,  Ln  =0.5  cm,  Mq  =  2.0  x  10~7  gm,  V^q 

=  6.0  x  107  cm/sec,  Z  =  10,  and  u  =  28.  In  Fig.  1  we  plot  the  position  of 

the  shell  R  (cm)  versus  the  time  t  (nsec).  Curve  A  shows  the  position  of 

a  'free  streaming'  debris  shell  and  is  based  upon  R(t)  =  V^t.  Curve  B 

shows  the  position  of  a  'decelerating'  debris  shell  and  is  based  upon  (7). 

The  two  curves  are  identical  for  t  <  20  nsec.  The  maximum  expansion 

radius  in  this  case  is  R  =  2.8  cm.  In  Fig.  2  we  plot  the  effective 

max  2 

gravitational  acceleration  gg££  (cm/sec  )  versus  the  time  t  (nsec).  It  is 

based  upon  (6).  We  note  that  geff  increases  substantially  during  the 

expansion.  In  Fig.  3  we  plot  yn/2.  versus  the  time  t  (nsec).  Here,  = 
1/2  u  i  u 

(ge  ^/L  )  is  the  growth  rate  of  the  magnetized  ion  Rayleigh-Taylor 

instability.  We  show  the  regions  of  stability  and  instability  based  upon 

the  criterion  Yq/W^  <  1  and  Yq/2^  >  1,  respectively.  We  point  out  that 

structure  was  observed  to  develop  at  t  =45  nsec  in  the  laser  experiment 

which  is  in  the  unstable  regime  (Ripin,  1987).  Finally,  we  present  a  plot 

of  Rg,  Rc,  and  RQgs  (cm)  versus  laser  energy  W  (joules).  In  calculating 

Robs  we  have  taken  a  =  1.5.  For  W  <  8  J  we  do  not  expect  that  structure 

caused  by  the  unmagnetized  ion  Rayleigh-Taylor  instability  will  develop 

since  R  <  R„.  However,  structure  may  not  be  observed  if  W  <  20  J  because 

C  D 

it  may  not  have  time  to  develop. 

The  structure  criteria  suggested  in  this  paper  [(12)  and  (13))  can  be 

checked  with  a  series  of  laser  experiments.  The  experiments  should  be 

designed  to  determine  the  scaling  of  the  observed  structure  radius  with 

relevant  physical  parameters  (e.g.,  y,  Mq,  Z),  and  to  determine  if 

structuring  stops  at  sufficiently  small  laser  energies  (i.e.,  W  <  W  ).  In 

1/2  c 

the  former  case,  the  scaling  Rq^s  «  Mq  would  probably  be  the  easiest  to 

verify.  Different  target  materials  could  be  used  to  scale  with  Z/y. 

In  the  case  of  aluminum  we  estimate  Z/y  =  10/28,  while  for  carbon  we 

expect  Z/y  =  1/2.  Much  heavier  targets  (e.g.,  gold)  may  yield  values  of 

Z/y  <  10/28.  Thus,  there  may  be  a  sufficient  range  Z/y  to  determine  if 

Robs  sca^es  linearly  with  Z/y.  Of  course,  there  is  uncertainty  in 

determining  Z  which  could  hinder  this  analysis.  On  the  other  hand,  if  the 
1/2 

scaling  with  Mq  is  verified,  lending  credence  to  (9),  then  we  could  use 
the  experimental  measurement  of  RQgs  vs.  Z/y  as  means  to  estimate  Z  for 
different  target  materials.  With  regard  to  checking  if  structure  is  not 
observed  for  sufficiently  small  laser  energies,  it  is  recommended  that  the 


parameters  Bq  »  10  kG  and  Z/p  =  10/28  are  used,  and  to  reduce  the  laser 

energy  to  10  J.  It  could  then  be  checked  if  structure  is  observed  for  V  > 

V  and  not  observed  for  V  <  V  . 

c  c 

IV.  DISCUSSION 

Ve  have  proposed  a  set  of  criteria  for  the  onset  of  structure  caused 

by  the  unmagnetized  ion  Rayleigh-Taylor  instability  in  magnetically 

confined  plasma  expansions.  The  AMPTE  magnetotail  barium  release 
(Bernhardt  et  al.,  1987)  and  the  recent  NRL  laser  experiments  (Ripin  et 
al.,  1987)  are  excellent  examples  of  such  expansions.  In  both  cases, 
structure  was  observed  on  time  scales  fast  compared  to  the  ion  cyclotron 
period,  and  the  unstable  density  shell  exhibited  'clumping'  as  opposed  to 
'rippling'.  Both  these  features  are  consistent  with  linear  theory  (Hassam 
and  Huba,  1987a, b)  and  2D  MHD  simulations  (Huba  et  al.,  1987). 
Furthermore,  the  onset  criteria  proposed  is  also  consistent  with  the  AMPTE 
and  laser  observations.  Aside  from  predicting  R£  (the  critical  radius  at 
which  structure  can  develop),  we  also  find  that  Rc  is  independent  of  the 
magnetic  field.  This  latter  point  has  recently  been  observed  in  the  NRL 
laser  experiment.  Finally,  we  have  recommended  a  series  of  experiments 
which  will  help  confirm  (or  refute)  these  onset  criteria.  Of  particular 
interest  is  the  scaling  of  R£  with  Mq,  and  the  onset  of  structure  for  W  > 

V  . 

c 

Finally,  we  comment  that  a  critical  parameter  in  determining  RQbs  is 
a.  At  this  point,  the  values  assumed  for  a  have  guided  by  observational 
data.  Although  Rc  is  independent  of  Bq,  it  may  be  that  a  is  dependent  on 
Bq  (e.g.,  through  the  growth  rate)  so  that  R0^s  could  be  dependent  on  Bq. 
For  example,  we  find  that  the  linear  growth  rate  of  the  unmagnetized 
Rayleigh-Taylor  instability  is  sufficiently  fast  to  explain  the  observed 
structure  position  for  both  AMPTE  and  the  Bq  =  10  kG  laser  experiment. 
However,  for  the  Bq  =  1  kG  laser  experiment,  the  growth  rate  of  the 

instabilty  (at  the  observed  wavelength)  is  too  slow  to  account  for  the 
structure  at  the  observed  position.  We  note  that  the  growth  rate  is 
proportional  to  the  wavenumber  (k)  or  mode  number  (m)  so  that  rapid  growth 
of  short  wavelength  modes  is  anticipated.  We  speculate  that  there  may  be 
a  rapid  nonlinear  cascade  of  energy  from  the  short  wavelength  modes  to  the 
long  wavelength  modes  that  are  observed.  We  are  currently  investigating 
this  hypothesis. 
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Plot  of  R  (cm)  versus  time  t  (nsec).  The  parameters  considered 
are  BQ  =  10  kG,  Mq  =  2.0  x  10 gm,  Ln  =  0.5  cm,  =  6.0  x  10^ 
cm/sec,  Z  =  10,  and  u  =  28.  Curve  A  is  for  a  'free  streaming' 
debris  shell,  while  curve  B  is  for  a  'decelerating'  debris  shell 
and  is  based  upon  (7). 
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Fig.  3  —  Plot  of  versus  time  t  (nsec).  The  parameters  are  the  same 

as  in  Fig.  1.  The  stable  and  unstable  regions  are  defined  by 
Tq/Bj  <  1  and  Tq/Bj  >  1,  respectively. 
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735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
OICY  ATTN  P.  FISCHER 
OICY  ATTN  W.F.  CREVIER 
OICY  ATTN  STEVEN  L.  GUTSCHE 
OICY  ATTN  R.  BOGUSCH 
OICY  ATTN  R.  HENDRICK 
OICY  ATTN  RALPH  KILB 
OICY  ATTN  DAVE  SOWLE 
OICY  ATTN  F.  FAJEN 
OICY  ATTN  M.  SCHEIBE 
OICY  ATTN  CONRAD  L.  LONGMIRF. 

OICY  ATTN  B.  WHITE 
OICY  ATTN  R.  STAGAT 

MISSION  RESEARCH  CORP. 

1720  RANDOLPH  ROAD,  S.E. 

ALBUQUERQUE,  NM  87106 
OICY  R.  STELLINGWERF 
OICY  M.  ALME 
OICY  L.  WRIGHT 


MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 

OICY  ATTN  W.  HALL 
OICY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 

I  2  3*40  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

OICY  ATTN  E.C.  FIELD,  JR. 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

OICY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

II  ARROW  DRIVE 
WOBURN,  MA  01801 

OICY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

OICY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  9H610 
ATTN  A.  THOMSON 

RAD  ASSOCIATES 

P.O.  BOX  9695 

MARINA  DEL  REY,  CA  90291 

OICY  ATTN  FORREST  GILMORE 

OICY  ATTN  WILLIAM  B.  WRIGHT,  JR. 

OICY  ATTN  WILLIAM  J.  KARZAS 

OICY  ATTIV  H.  ORY 

OICY  ATTlf  C.  MACDONALD 

OICY  ATTN  BRIAN  LAMB 

OICY  ATTN  MORGAN  GROVER 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY  ,  MA  01776 

OICY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
3  30  WEST  «2n<3  STREET 
NEW  YORK,  NY  10036 

OICY  ATTN  VINCE  TRAPANI 


SCIENCE  APPLICATIONS 

INTERNATIONAL  INCORPORATED 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 

0 1  C  Y  ATTN  LEWIS  M.  LINSON 
01CY  ATTN  DANIEL  A.  HAMLIN 
01  CY  ATTN  E.  FRIEMAN 
0 1 C Y  ATTN  E.A.  STRAKER 
01 C  Y  ATTN  CURTIS  A,  SMITH 

SCIENCE  APPLICATIONS 

INTERNATIONAL  CORPORATION 
1710  GOODRIDGE  Dft. 

MCLEAN,  VA  22102 
0 1 C  Y  J.  COCKAYNE 
01  C  Y  E.  HYMAN 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  9^025 


01  CY 

ATTN 

J.  CASPER 

0 1  C  Y 

ATTN 

DONALD  NEILSON 

01  CY 

ATTN 

ALAN  BURNS 

0  1  C  Y 

ATTN 

G.  SMITH 

01  CY 

ATTN 

R.  TSUNODA 

0 1  C  Y 

ATTN 

DAVID  A,  JOHNSON 

01  CY 

ATTN 

WALTER  G.  CHESNUT 

0  1  C  Y 

ATTN 

CHARLES  L.  RINO 

01CY 

ATTN 

WALTER  JAYE 

0 1  C  Y 

ATTN 

J,  VICKREY 

0 1  C  Y 

ATTN 

RAY  L.  LEADABRAND 

0  1  C  Y 

ATTN 

G,  CARPENTER 

01  CY 

ATTN 

G.  PRICE 

0  1  C  Y 

ATTN 

R.  LIVINGSTON 

01  CY 

ATTN 

V,  GONZALES 

0  1  C  Y 

ATTN 

D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 

75  WIGGINS  AVENUE 

BEDFORD 

,  MA  01730 

0  1  C  Y 

ATTN 

W.P.  BOQUIST 

TRW  DEFENSE  & 

SPACE  SYS  GROUP 

ONE  SPACE  PARK 

REDONDO 

BEACH 

,  CA  90278 

0  1  C  Y 

ATTN 

R.  K.  PLEBUCH 

01  CY 

ATTN 

S,  ALTSCHULER 

0  1  C  Y 

ATTN 

D.  DEE 

01  CY 

ATTN 

D/  STOCKWELL 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MA  01803 
01 CY  ATTN  W.  REIDY 
01CY  ATTN  J.  CARPENTER 
01 CY  ATTN  C.  HUMPHREY 

UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 

0 1 C  Y  ATTN:  N.  ZABUSKY 


i 


ft 


